GAUSSIAN FIELDS AND GAUSSIAN SHEETS WITH 
GENERALIZED CAUCHY COVARIANCE STRUCTURE 
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Abstract. Two types of Gaussian processes, namely the Gaussian field with 
generalized Cauchy covariance (GFGCC) and the Gaussian sheet with general- 
ized Cauchy covariance (GSGCC) are considered. Some of the basic properties 
and the asymptotic properties of the spectral densities of these random fields 
(-H ^ are studied. The associated self-similar random fields obtained by applying 

-4-^ , the Lamperti transformation to GFGCC and GSGCC are studied. 

1. Introduction 

Generalization of some well-known stochastic processes indexed by a single pa- 
rameter to processes indexed by two parameters has attracted considerable interest 
recently. For example, Houdre and Villa [25 generalized fractional Brownian mo- 
p^ I tion parametrized by a single Hurst index to the bifractional Brownian motion 

f^ . characterized by two indices. Another example is provided by the multidimen- 

^D ' sional stationary Gaussian fields with generalized Cauchy covariance indexed by 

[~^ . two parameters introduced by Gneiting and Schlather [22] . These processes can be 

^^ I regarded as extension of the Gaussian processes with Cauchy covariance used in 

OO . geostatistics. For simplicity, we call such processes the Gaussian field with general- 

ized Cauchy covariance (GFGCC). Here we would like to point out that one should 
not confuse such a process with a stable process with Cauchy marginals. 

In general, processes parametrized by two indices can provide more flexibility 
in their applications in modeling physical phenomena. In particular, the GFGCC 
C^ ' model has an additional nice and useful property as it allows separate charac- 

terization of fractal dimension and long range dependence (LRD) by two different 
indices. This is in contrast to models based on fractional Brownian motion (or frac- 
tional Brownian noise) which use a single index to characterize these two properties. 
Models based on a stochastic process or field parametrized by a single index seem 
inadequate. A detailed analysis on network traffic carried by Park et al [H] shows 
that fractional Brownian motion/fractional Gaussian noise model is inadequate for 
description of network traffic for all scales since at very small time scales the traf- 
fic fluctuations are no longer statistically self-similar. The need to replace global 
scaling by local scaling is essential for processes such as multifractional Brownian 
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motion introduced independently by Peltier and Vehel [45] and Benassi, Jaffard 
and Roux [7]. 

The main aim of this paper is to study GFGCC and its anisotropic counterpart, 
which we call Gaussian sheet with generalized Cauchy covariance (GSGCC). In 
view of the fact that GFGCC is widely used in geostatistics and other applications 
[nilHlllHlEIlEOlllOllllllin], it will be useful to consider its properties in more detail. 
In these existing applications, usually only the covariance structure of GFGCC are 
used, and the sample path properties of GFGCC are rarely mentioned. However, a 
better understanding of the sample properties of GFGCC and GSGCC will render 
more versatility and flexibility to their applications. Our approach to this subject 
is mainly from a physical viewpoint. Basic sample properties such as the long 
range dependence and the local self-similarity properties of GFGCC and GSGCC 
are investigated. A simpler method is used to derive the asymptotic properties 
of the spectral densities of GFGCC and GSGCC. By generalizing the Lamperti 
transformation to n-dimensional processes, new types of random field and random 
sheet with global self-similar property associated with GFGCC and GSGCC are 
obtained. Properties of these random field and random sheet are also studied. 

2. Isotropic Gaussian Field with Generalized Cauchy Covariance 

In this section we consider GFGCC, which is a multidimensional isotropic Gauss- 
ian random field in n-dimensional Euclidean space. We first introduce some no- 
tations and state some basic definitions and properties of GFGCC. Denote by N, 
Z, R and K+ the sets of positive integers, integers, real numbers and positive real 
numbers respectively. Let t = (ti, . . . , t„) and s = (si, . . . , s„) be two vectors in 
M", and ||t|| = \/X]"=i ^h '*'J e N be its Euclidean norm. By t — > 0+ and i — > oo, 
we mean U -^ 0+ and ti ^ oo respectively for alH = 1, . . . , n. 

Definition 2.1. A random field Xa,f3{t) on M" is called a Gaussian field with 
generalized Cauchy covariance (or GFGCC) if it is a stationary Gaussian field with 
mean zero and covariance given by 

(2.1) C^^^ir) = (X,,^(t + T)X„,^(t)) ^{1 + ||r|r)-^ , 
where a G (0, 2] and /? > 0. 

Note that ()2.1|) has the same functional form as the characteristic function of the 
generalized multivariate Linnik distribution first studied by Anderson [2]. Ca^isir) 
is positive-definite for the above ranges of a and /3, and it is completely monotone 
for < a < 1, /3 > 0. Xa.p{t) becomes the Gaussian field with usual Cauchy 
covariance when a = 2, /3 = 1. 

Recall that a random field X{t) is iJ-self-similar (Hss) if X{ct) =d c^X{t), 
where =d denotes equality in the sense of finite-dimensional distributions of X. 
Self-similar property requires scale invariance to hold for all scales. This is rather 
too restrictive for many applications. We also know from Samorodnitsky and Taqqu 
[47] that a stationary Gaussian random field such as Xa,p{t) can not be a self-similar 
field. However, Xa^pit) satisfies a weaker self-similar property known as local self- 
similarity considered by Kent and Wood [28]. 

Definition 2.2. Let a e (0,2]. A centered stationary Gaussian field is locally 
self-similar (Iss) of order a/2 if for ||t|| -^ 0+, its covariance C{t) satisfies 

(2.2) Cir)=A-B\\rr[l + 0{\\r\f)] 
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for some positive constants A, B and S. 
Since as ||t|| -^ 0+, 

(2.3) Co..Ar) = l-(3\\r\ni + \\r\n, 

the GFGCC Xaj3{t) is a/2 - Iss, with A = 1, B ^ f3 and 6 = a. Adler [1 called the 
class of Gaussian fields which satisfy (|2.2p the indexed-a fields. These processes are 
also known as the Adler processes according to some authors, for example Lang and 
Roueff |34| . They form a very rich class of Gaussian random fields, which include 
the centered Gaussian field Sct.^(t) with powered exponential covariance 

(2.4) (S„,Mt + T)S„^^(t))=e-^ll^ll°, 

which have the same functional form as the characteristic function of the multi- 
variate symmetric stable distribution as given in Kotz, Kozubowski and Podgorski 
[31] , and Garoni and Frankel [TTj . Instead of using (|2.2p to characterize local self- 
similarity, one can also use the definition of locally asymptotically self-similar (lass) 
property first introduced by Benassi, Jaffard and Roux [7j for multifractional Brow- 
nian motion Bjju\ (t) , which is a generalization of fractional Brownian motion with 
the Hurst index replaced by the Hurst function H{t), < H{t) < 1. It can be 
shown that under some regularity conditions on H{t), the multifractional Brown- 
ian motion -Bff(t) (t) is lass. This property can be adapted to GFGCC if we take 
H{t) as constant with its value in (0,1). 

Definition 2.3. A stochastic process X{t) is lass at a point to with order k if 

2.5) hm <^ \ =d Tt„{u), 

and Ttg (u) is nontrivial. Here the convergence and equality are in the sense of finite 
dimensional distributions, and Ttg{u) is called the tangent field of X{t) at the point 
to- 

Proposition 2.4. GFGCC is a lass random field of order a/2; and its tangent 
field is Levy fractional Brownian field of index a/2. 

Proof. By using (|2.3p . the covariance of the increment field ArXa.pit) :— Xa.pit + 
t) — Xa,i3{t), for p, CT ^ 0+ is given by 

{A,x^At)A.x^At)) = p (iipir + ikir - \\p <^\n + o (hph^-, ikn^-) . 

Let p — eu and a — ev, then 

,. / A.suXa^i3{t) Ag^Xg^pit ) ' 

e ™+ \ e"/2 £"/2 

= lim {/3(||M|r + ||i'r-||u-u||") + 0(e")} 

=2/3(B,/2(u)B„/2(i;)), 

where -Bq/2(u) is the index-a/2 Levy fractional Brownian field with zero mean and 
covariance given by 

(2.6) {BgA^)Bg/2{v))^\(\\ur + \\vr-\\u^vr), u,veK^. 

Thus up to a multiplicative constant \/2p^ the tangent field of GFGCC at any point 
io € M" is the Levy fractional Brownian field indexed by a/2. D 
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The tangent field of Xa.jsit) at the point to reflects the loeal structure of the 
random field at to. In other words, GFGCC behaves locally like a Levy fractional 
Brownian field. This provides an example to the general results on tangent fields 
considered by Falconer [15j . 

The fractal dimension of the graph of a random field X{t) depends on the local 
property of the random field. The local irregularities of the graph are measured 
by the parameter a, which can be regarded as the fractal index of the random 
field. Thus the behavior of the covariance function at the origin to a great extent 
determines the roughness of the random field. The results on the fractal dimension 
of an Iss field are treated by Adler [1], Kent and Wood J28], Davies and Hall 12j. 
A nice property of a/2-lss fields is that their fractal dimension is determined by a. 

Definition 2.5. Let X{t) be a stationary Gaussian field and let (j'^{t) — (^ArX{t)'^'j 
be the variance of the increment process ArX{t) := X{t + r) — X{t). If there exists 
a e (0,2] satisfying 

(2.7) a=sup{(^ : cr2(T) =o(||T|n as ||r|| -> O} 

= inf{^ : \\T\\'^ = o{a^{T)) as ||r|| ^ 0} , 

then a/2 is called the fractal index of the random field X{t). Equivalently, a/2 is 
the local Holder index of the random field. 



Clearly, condition (|2.7p is fulfilled by random field which satisfies (12.21) such as 
the GFGCC Xa,p{t). Adler had shown that a/2 is the upper bound of the indices 
for which, with probability one, the graph oi X{t) satisfies a global regularity of the 
same order. Thus, a characterizes the roughness of the sample path. The fractal 
dimension of GFGCC can be obtained by using the following result for the fractal 
(Hausdorff) dimension of an Iss field as given in Adler 1 , Chapter 8. 

Proposition 2.6. The fractal dimension D of the graph of a locally self-similar 
field X{t),t € R", of fractal index a/2, over a hyperrectangle C — nr=i[°'*' ^*]' *'^ 
given by 

(2.8) D^n+1-^. 

The estimation of a for Iss field has been studied extensively, see for example 
Wood and Chan [5T], Istas and Lang [2^, Kent and Wood 28J, Lang and Roueff 
|34j . The parameter /3 is also known as topothesy in the studies of roughness of 
surfaces by Wood and Chan 51j, and Davies and Hall [12,. The topothesy of a cross 
section provides a measure of the roughness which is scale-dependent in contrast 
to fractal dimension which is scale-invariant. 

The Gaussian random field Xa,i3{t) can have SRD (short range dependence) or 
LRD (long range dependence), depending on the values of the parameters a and (3. 
For this purpose we make use of the following definition which is a generalization 
of the one-dimensional case considered by Flandrin et al ^16J , Lim and Muniandy 

m- 

Definition 2.7. A stationary centered Gaussian field with covariance C{t) is said 
to be a long range dependent process if 

(2.9) / |C(r)|d"r = oo. 
Otherwise it is short range dependent. 
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Proposition 2.8. The GFGCC Xa^pit) is a long range dependent random field if 
and only if < al3 < n. 

Proof. In order to obtain the condition for the Gaussian random field with co- 
variance (j2.f p to be LRD, we make use of the foUowing integral identity given in 
Gradshteyn and Ryzhik i23j, 3.251, no. 11: 



/•OO 1 

/ x^'-^il + xPydx^-B 
Jo P 



P P 



where p > 0, < /i < pw and B{x,y) — r{x)T{y)/r{x + y) is the beta function. 
Using polar coordinates, one gets 

(2.10) / \C^Ar)\d-T=l (l + ||r|r)-^d"r 



+ 
27rt 



2«r(§) 







CX). 



For large r. 

Therefore, the integral (|2.10p is divergent for all /3 > 0, < a/3 < n. For /3 > and 
aP > n, we have 

Therefore the condition for Xa.pit) to be a Gaussian field with LRD is < a(3 < 

n. D 

The discussion above shows that it is possible to characterize the fractal di- 
mension D and the LRD property separately. If the covariance is re-expressed as 
(1 -I- ll'T'll")"''' , which behaves like ||t||~''' in the large-||T|| limit, then ^^^^/^(i) 
is LRD if and only if < 7 < n. Thus a (0 < a < 2) and 7 (7 > 0) respec- 
tively provide separate characterization of fractal dimension and LRD/SRD. The 
separate characterization of the fractal dimension (local property) and LRD (global 
property) for GFGCC appears to offer a more natural and flexible model than that 
based on a single parameter such as in Levy fractional Brownian field. We note 
that this feature of separate characterization of local self-similarity (hence fractal 
dimension) and long range dependence is present in any stationary Gaussian field 
with covariance C{t) satisfying the asymptotic behaviors C{t) ^ A ~ _B||t||" as 
||t|| -^ 0+, and C{t) -^ \\T\\-'^' as ||r|| -^ 00, with a G (0,2], 7 > 0. Similarly, one 
can also have a Gaussian stationary process which has separate parametrization of 
fractal dimension and short range dependence [39]. The ability to have separate 
characterization of fractal dimension and Hurst effect is a desirable property in the 
modeling of physical and geological phenomena. 

3. Asymptotic Properties of Spectral Density of GFGCC 

In this section, we consider the spectral density of Xa,i3{t) and its asymptotic 
properties. Though the covariance of GFGCC is given by a relatively simple ex- 
pression, the analytic simplicity of the covariance function is not inherited by the 
corresponding spectral density. This is similar to the case of the stationary Gauss- 
ian field Sa,/3(t) with powered exponential covariance (|2.4p which have simple form. 
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but its spectral density in general does not have closed analytic expression. In the 
case of GFGCC in R, a detailed study of spectral densities (in terms of probability 
distributions correspond to the characteristic functions of generalized Linnik dis- 
tributions) have been carried out by Kotz et al. 30J for < a < 2, /3 = 1, n = 1; 
by Ostrovskii [42] for < a < 2, /3 = 1, n €E N; and by Erdogan and Ostrovskii 
[M] for 0<a<2,/3>0,n = l. They employed the contour integration repre- 
sentations and series expansions of the generalized Linnik distributions. However 
the techniques used in these works are less accessible to practitioners. In this sec- 
tion, we derive the asymptotic properties of the spectral densities of GFGCC for 
< a < 2 and (3 > 0, which can be regarded as an extension to the results on gen- 
eralized multivariate Linnik distributions. The techniques used in our derivations 
are mathematically more tractable. 

Recall that the spectral density S{uj) of a stationary field X{t) is defined as the 
Fourier transform of its covariance function C{t) ~ {X{t)X{0)): 

^(")=(2^/^"^"'^(')'^"'' 

if the integral is convergent. If the integral does not converge, we consider C{t) 
as a generalized function and define S{uj) as the Fourier transform of C{t) in the 
Schwartz space of test functions [18 . Namely, for any test function iP{lu) in the 
Schwartz class of M" , we require 

where 

^^(*)=(2^/^"^^'^(")'^""- 
Alternatively, the spectral density can also be defined to be the function satisfying 

For GFGCC, since 

(3.1) j,(,)^^_cos(z----j 

as z — > oo ([3], page 209), we find that when a[3 > ^^^-j^, it's spectral density is 



2-71 

\UJ\\ 2 



J^(ll^llt) 



''^ (27r)"i (l + ||t|l")/3 (27r)2 ig (1 + i")/3 

where J,y(z) is the Bessel function. When a = 2, using no. 4 of 6.565 in [53], we 
have the explicit formula 

(3-3) ^..M-) - ^lJMS?(^^t-^(ll-ll)' 
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ii (3 > {n — l)/4. Here K^{z) is the modified Bessel function. The formula no. 7 of 
6.576 in [25 shows that if /3 G (0,7i), then 

V2t+''-i7rtr(/3) ^ '^^" "7 

Therefore (|3.3p is stiU the spectral density when /3 e (0, (n — l)/4]. For general 
a < 2, no explicit formula such as (|3.3p can be found for 5'q_^(w). When n = 1, 
the formula p.2p gives the spectral density of Xa^pit) for all values of a e (0,2] 
and /3 > 0. For n > 2, we would also like to find a formula for the spectral density 
that is valid for all a € (0, 2) and /3 > 0. For this purpose, it would be beneficial 
to investigate the case n—1 first. When n = 1, we can rewrite p.2p as 

Let 

and consider the region S^ in the complex plane defined by 

S)r = |z e C : \z\< r, Rez > 0, Imz > o|. 

When a G (0, 2), the function / is an analytic function on the domain S^. There- 
fore, by Cauchy integral formula, 

(3.4) ^ fiOdC = 0. 

Notice that the boundary of T)r, 31) r, consists of three components: the line seg- 
ment Ir^i along the real axis from to r, the arc Cr of the circle \z\ = r from r to 
ir, and the line segment 2^,2 along the imaginary axis from ir to 0. On the arc Cr, 
if r > 1, then 

-|w|ImC 

1/(01 < 



Therefore 



and (|3.4|) implies that 



lim / /(C)dC = 0, 



r^oo 



hm / /(C)dC = - lim / fiOdC 
This gives us 
(3.5) Sa,p{uj) =-Re / , ,g = Im / -^du. 
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For n > 2, we can derive a formula similar to ()3.5p . Recall that the Hankel's 
fmiction of the first kind He, (z) is defined as 

Hi^Hz) = Mz)+iN,{z), 

where N^(z) is the modified Bessel function of the second kind or called the Neu- 
mann function. Using Hankel's function, we can rewrite p.2p as 

II 11^ .^iJ«,(||^||i) 

For z -^ cx), we have (|13], no. 3 of 8.451) 



Hl'>{z) 



2 r / Tiv t:^ 

Therefore, we can show as in the n = 1 case that 



<S'a,/3(w) = - 



{27r)^ 



H^ .^ ill:2.(z||a;h) 



■Im 






i-(iu) 2 (iu. 



Using the formula ([23 , no.l of 8.407) 



^.W = ye^ijW(zz), 



we have finally 

(3.6) Sa.pit^) 



\io\\^ 



2 2 TT 2 



Im 



if„-2(||a)||M) 



° [l + e^M"] 



u2 (iu. 



This formula agrees with the formula for multivariate Linnik distribution proved 
in [42] for a e (0, 2), /3 = 1 and n e N. Notice that the right hand side of ((a!6|) is 
well-defined for all a,/? > 0. Using the formula (no. 2 of 6.521 in [23]), 



xK,j{ax)J^(hx)dx 



f 



a-'ia^ + b^)' 



V > -1, 



we have 



-Im 



2 2 TT 2 Jo 



f l + e^-u") 



=(27r)2||i||T- / J„-2(^||t||)tj^ 
-'0 ^ 



LO 2 



^i^ Il±2 
2 2 TT 2 



u^ du 



Im 



W"w 



-Rrn-2 (u;m) 



tU^ du } duj 



M^lm ^"^ 



u: 



-Im 



1 + e 2 yC 



^ -^0 (l + e'¥^U")^(M2 + ||f||2) 



-J / LLjKn-2 {uJU]J-n-2 { Uj\\t\\) dUJ du 

-du 
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= : / ^ du. 

When a E (0, 2), residue calculus implies that this last integral is equal to 

u _ I 

(M-i||t||)n+e^M"j ^^ + 11^11 I 

This shows that p.6p is indeed the spectral density of GFGCC for all a G (0, 2) 
and /3 > 0. We would also like to remark that although the formula (|3.6p is derived 
under the assumption n > 2, but since 

therefore when n = 1, the formula p.6p reduces to the formula (|3.5p . We summarize 
the result as follows. 

Proposition 3.1. If a E (0,2) and /? > 0, the spectral density of the GFGCC 
Xa,0{t) is given by 

IIljII^ r°° Kr,-2 (\\lu\\u) 

(3.7) S^^u;) = -Jj-^-^Im / ^ -^u^du. 

2 2 TT 2 Jo ^i + e^^a^ 



j oo 



If a ~ 2 and /3 > 0, i/ie spectral density of the GFGCC X2.p{t) is given by 

To find the high frequency behavior of the spectral density, we first consider the 
case where a = 2. Using the fact that ([23], no. 6 of 8.451) as 2 — > 00, 

2z^ ^(2z).-j-!r(^-j + i)- 
This implies that as ||aj|| — > cx), 

^'■'^ ^^■'^^"^ 2-^+%-^r(/3)' ^^0 (2|MIV-.!r(H±i_,„;3)' 

with leading term 

ll„,ll/3-^ „ „ 

5'2,/3(w) 



2^+/37rT-r(/3) 



For general a S (0, 2), to find the high frequency behavior of Sa.p{io)^ we make 
use of eq. p.7|) . Making a change of variable and using 



(3.10) ^ — -, = E ^^^^^^^ + o(ii.ir"(-^)). 



i + ^^frl -° 
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as 


a; 


-^ 


oo, 


we find that 


Dq- 


A^) = 


= - 




11/ Ji — ^ 


2' 


1-2 ,1 + 2 J-™ 
2 TT 2 






= 


= - 




IL ill —n 
\\UJ\\ 




2" 


1-2 „ + 2 -p/ 
2 TT 2 i 1, 



wt 



K n-2 (U) — nr'^" 

^ ' 



22 W 



[B + j)e 2 - — - — -u^du 



+ 0(||cj||-°(™+i)-") as||tj||^oo. 

Using tlie formula ([23 , no. 16 of 6.561) 

j^ x^K,ix)dx = 2^-ir (^if^) r (^^f^) , Re (/i + 1 - IH) > 0, 

and tiie definition of asymptotic expansion ([3], page 611), we conclude that as 
||a;|| -^ oo, Sa,p(oj) behaves asymptotically as 

(3.11) 

Notice that there is a drastic change of high frequency limit of Sa,f3{io) when a < 2 
and a = 2. In fact, naively putting a = 2 in p. lip give identically zero terms. 
This is a hint that as ||a;|| -^ cx), 5*2, /5(w) does not have polynomial decay, instead 
it decays exponentially as is verified by (|3.9p . We summarize the results as follows: 



Proposition 3.2. If a € (0, 2) and f3 > 0, the high frequency limit of the spectral 
density Sa,/3ioj) is given by the following asymptotic series 

(3.12) 

5.,M-)-;iT2^E^^-^r(/3 + ,)r(^^jr^-^jsm^lMI - ". 

If a — 2 and /3 > 0, the high frequency limit of the spectral density S2,fi{oj) is given 
by the following asymptotic series 






2 



2^+%^r(/3) ^ (211^11)^" j!r (ii±i - J- - /3) ■ 

When a e (0,2),/5 = l,n G N and a € (0,2),^ > 0,n = 1, ((XT2l) agrees with 
the results given in [42^ and 14] respectively. In particular, we observe that when 
a e (0, 2), the high frequency behavior of the spectral density of GFGCC is 

(3.13) 

ScA^) - ^r ^— r ^— sin— ||a>|r"-" ^ 0+, llc^ll ^ oo, 



TT 2 



2 



which is independent of /3. Kent and Wood [55] have shown that if a random field 
has spectral density satisfying (|3.13p , then its covariance satisfies (|2.2p with locally 
self-similar property. However, the converse is not true. 
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In this connection we remark that for the Gaussian stationary field SQ,j^(t) with 
powered exponential covariance which is Iss with 

(3.14) 

i^c^At + T)S.,Mi)) - e-'^ll-ll" = 1 - ;3|lr|r [1 + 0(|lr|r )] , |lr|l ^ 0+, 



its small ||r|| behavior has a similar form as that of GFGCC (|2.3p . Thus it is not 
surprising that this two random fields have the same tail behavior for their spectral 
densities at high frequencies as given by (jS.lSp . The detailed calculation carried out 
by Garoni and Frankel [T7] for the probability distribution of the multivariate Levy 
stable distribution with characteristic function given by (|3.14p confirms this. We 
also note that (I3.12p can be used to verify that the tangent field at any point to has 
spectral density which varies as |jijj|j~"~" for ||a;|| -^ cx). If we let a = 2H, then the 
tangent field is just the Levy fractional Brownian field in M". Such a relationship 
can be viewed as a consequence of the Tauberian-Abelian theorem (see e.g. [29 ). 
For the low frequency behavior of the spectral density Sa^/j^uj), we first consider 
the case a = 2. Using 8.485, 8.445, 8.446 of [23], we find that if j/ ^ Z, 

(3.15) K^z) K^.[z) 23i,(,,) jl. ^.,r(, + 1 _ ,) 1. ,,r(, + 1 + ,) 

whereas when v — ±m, where ttt, is a nonnegative integer, 
1 "^ {-ly {m - j - ly. /z\2j-™ 



(3.16) K,{z) =- J2 



2 -^ j! V2 



Here ip{z) — r'(z)/r(z) is the logarithm derivative of the Gamma function. There- 
fore from (13.81). we find that: 



r(/3- 



• if /3 > n/2, then as w 


— s 


■o+. 


(3.17) 




S2A^) ^ 


• if /3 = n/2, then as tj 


-) 


■o+. 


(3.18) S2,i3{^) - 


^¥- 


1 

^-i7rtr(f) 

1 



2«7r-r(/3) 



{-ln||cj||+ln2 + V-'(l)} 



„ , , {-ln||cj|| +ln2-7|. 

where 7 is the Euler constant. 
• if /? < n/2, then as \\lu\\ -^ 0+, 
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Figure 1: The spectral density 5'a,/3(||aj||) as a function of ||lj|| when n = 3 and a/3 = 1.5. 



In fact, by considering the cases /3 — ^ G Z and /3 — ^ 7^ Z separately and 
substituting the series (|3.15p and (|3.16p into (|3.8p . we can express the spectral 
density 52, /^(w) in terms of convergent power series in ||ci;||. 

For general a G (0,2), the low frequency behavior of Sa.fiii^) depends on the 
arithmetic nature of a and /3. The method we are going to employ does not allow 
the derivation of the whole asymptotic series as obtained by Kotz et al [3U] , Erdogan 
and Ostrovskii [42l [14] . We will only derive the leading behavior of the spectral 
density Sa,i3{(^), which only depends on the algebraic conditions af3 > n, a(3 = n 
or ck/S < n. These conditions are less stringent than the arithmetic conditions 
considered in [30l|42l[T4]. However, the simpler method employed here provides 



GAUSSIAN PROCESSES WITH CAUCHY COVARIANCE 



13 











; \ \ 


a=0.4 


0.08 


' \ ' 


a=0.8 - 






a=1.2 


0.07 


I \ ^ 


a=1.6 . 




\ \ 


a=2 


0.06 


* \ \ 


- 




* \ '' 


f 005 


' \^ 


'~ix 


\ \ \ 


05° 0.04 


\ \\ 




\ \ N 




s \ ■> 


0.03 


^ \ ^ 














0.02 










I ■■', ~" ~, ^"^-^^Cj" ~ 






0.01 


"\ ''•■., ~'~. ~ ^---C_;:_~ ~ 




^-_ __ ~'~'~y^Z'?^~r^^^~=^- 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 





xlO 


































1 V 
























a=0.4 


0.9 




\ y. 
























a=0.8 


0.8 




1 V 
























a=1. 


2 






0.7 




























01=2 








0.6 






























3" 
==.0.5 

0.4 




H 


























0.3 




\^ 


^h 
























0.2 

0.1 

n 


-^^ 


~ 7 r 






s^ 


^v 


^ 


_^ 













2 3 4 5 6 



7 8 

IIHI 



10 11 12 



Figure 2: The spectral density S'a,/3(||(x;||) as a function of ||cij|| when n — 3 and q/3 — 3. 



the necessary Saf3{uj), \\uj\\ — )■ 0^ asymptotic behaviors which are sufficient for most 
practical purposes. 

When a/3 > n, since ([2?;, Eq. 8.402) 



(3.20) '^^W-Ff(^7TT)+^(^^"') 

we find from p.2p that as ||tj|| -^ 0+, 
(3.21) 

1 /-OO JL71 — 1 

Sa.jsil^) ^ : — TT^-r-T- I — —rrrdt 



as z ^ 0, 



1 r(^)r(/3-^) 



2"-i7rtr(f) 7o {l+t°'f 2"-i7rtr(f) aT{p) 
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Figure 3: The spectral density Sq,,/3(||(.<j||) as a function of \\u}\\ when n = 3 and a/3 = 4.5. 
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], no. 11 of 3.251). When a (3 = n, we make a change of variable on (|3.2p to get 



t2 



'^^ ' (27r)t io ^^ ^(||t^|h + i")'3 
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In view of the leading behavior of the Bessel function J^{z) as z ^ p.20p . we 
write Sa,/3{uj) as the sum of two terms S*^ ^(w) and S^ o{uj) where 

Jn-2 (t)—r, n Ta-dt, 



(27r)t 7i -^^ ' {\\oj\\°'+t'^)0 ' 
and 

As ||aj|| -^ 0+, p.ip and (|3.20p show that S^ p{llj) has a finite hniit. By Lebesgue's 
dominated convergence theorem, the hmit is given by 5'^ ^(O)- Namely 

1 /"^ / t^ \ 1 f°° 

J-n-2{t) ^^2 t-^di+- -^ / Jri-2{t)t~^dt 



This expression can be evaluated using regularization method. More precisely, using 
Lebesgue's dominated convergence theorem again, we find that 

SlJO)= lim < 7-^ / ( Ji^W-^ri^ ]t^^^'dt + —^ f J:^{t)t-^+'dt 



0+ (2^)^7o V ~ 2^r(t)/ (2^) 



2y 

1 



= lim l—^ f J.~2(t)t-^+'dt , 1 , , / t-^+'dt 

e^o+\{27:)- J„ —^' 2»-i7rtr(f)yo 

The formula no. 14 of 6.561 in [53] then gives 

{2^r (-) 111 

2«^tr('2^)"2"-%tr(f)7j 



lim ; — jT — -. — 7- — 

6^0+ 2"-i7rtr(f) £ 



{(l+eln2)(l + £^(I))(l + £v^g))-l} 



^ — ^^ (2 In 2 + VXl) + V' f-) I 

2«7rtr(f) I ^^ ' ^ \2)l 



For the term 5^ rt('^)i we make a change of variable m = t" or equivalently t — 1 
to get 

"^M'^^ 2"-%tnr(f)yo (||^lh+w)/5- 
We split S"^ /3('^) again into a sum of two terms S^ Juj) and S** n{'-^)i where 

(I) ||tj||" 2"-i7r*r(f) llt^ll ^" " ^' 



2»-i7r-nrni 
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and 

/3 rj^uf^-^-{i + uf-^ 



'2"-i7rtnr(f) io {l + uY 

When \\lj\\ -f 0+, 

The integral is a convergent integral with value given by ([23], 3.219 page 316) 

r-oo ^/3-l _ fl + u)/'-! 

. (ii,y' ''"=-w)-7. 

Putting everything together, we find that 

Therefore, as \\u}\\ —f 0+, 
(3.22) 

When a/3 < n, care has to be taken since p.2p is defined only for a/3 > ■^^-j-'-. 
Making a change of variable, one finds that 

(27r)2 7(, (||w||" + u")P 
It is easy to verify that for any a > and /3 > 0, when \\u!\\ -^ O"*", 

Sa di^) ~ —, — TTi — lim lim / ,,^ , —r-du. 

(27r)t e-*o+||.^|Hoio (||t^|h + u")'3 

Now using no.l of 6.621 in [53], we have 
lim / ^if-r urdu^ J„-2 (u)e-''"M^-"^ 



oJo (ll'^lh + '^") 







1 r{n~a/3) (n-aP aP-l n 1 



■ 2 



Fi 



2"^ V(£' + l)"""^r(f) V 2 2 '2'l + £ 



r2 



where 2-^1(0, 5; c; z) is the hypergeometric function 

^ , , N x-^ (a)i(&)i -Z-' / s , s , N r(a: + 7) 

p; (c)j j! r(a;) 
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Using the formula no. 2 of 9.131 in [23], we find that 

n-aP a(3-ln^ 1 ^^ _ T (f ) T (i) ^ fn - aP a(3 - 1 I ^ e^ 



2-ri 



-^1 — 7i — ' — 7i — ; tt; "T": — 9 — — 7^~:\ — 7 — ~: — ::v 2^1 



2'1 + eV p/'c»/3\ p / »+i-a/3 \ 'V 2 ' 2 '2'l + e2 
r(f)r(-i) ^/a/3^ n- a/3 + 1.3. 



i + e^y r(^)r(^ 

Since 



lim 2-Fi(a,5;c;z) = 1, 




therefore, 




/•- J„-2 (u)e-^«ut ^ r(n-a/3) 


o*-Q/3 V "^ / 

r(t) 


e^0+ ll<^II^O /o f||w||" + u")'3 r,i^ ^ /a/3\ J, /„+i_„/3\ " 



where we have used the formula T{2z) = 2^^ ^tt 2r(z)r(z + (1/2)). This gives for 
< a/3 < n, 

(3.23) 5„_^(,,)^iJL_ yy as||c.||^0+. 

It is easy to verify that by putting a = 2 in p.2ip . (|3.22p and (|3.23p . we get back 
(|3.17p . (|3.18p and p.l9p respectively. The low frequency behaviors of Sa^piio) are 
summarized below. 

Proposition 3.3. For all a £ (0, 2] and (3 > 0, the low frequency limit of the 
spectral density Sa^p{uj) is given by 

n T^ ( n — afi\ 
||^||Q/3-n i I — 2"^ I 

SaA^) ^rr^ / ^x : ifaP<n; 

2°'Ptt2 P ( M ) 

^..M-)-^;rzi4^r(^{ln^-^(V.(/5)+7)+ln2-i7+^^Q)}, «/ «/? = n; 

S^A^) - 2»-r.tP(.) ^"Ug) ' '^ "'^ " "■ 

We see that if < a/3 < n, Sa,p{(^) is divergent at the origin and if < a/3 < n, 
Sa,p{^) ^ ll^ll"''"" as ||a;|| -^ 0^. Note that the condition < a/J < n agrees with 
the LRD condition. In fact, it is a basic fact that for a stationary field with positive 
covariance, it is LRD if and only if its spectral density diverges at the origin. We 
would also like to point out that the low frequency limit of the spectral density of 
the Gaussian field with powered exponent covariance 'E.o,p{t) (|2.4p is finite as shown 
by Garoni and Frankel 17 , in agreement with the fact that SQ,^^(i) is SRD. 

We remark that by considering 

Sa.iiit-o) — "leading order term as Ht^H -^ 0^", 

we can find the next order term in the asymptotic expansion of Sa./siuj) when 
||ti;|| — > 0^ using the same methods we employed above. The results depend on 
more complicated conditions on a and f3. 
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We also briefly remark that the asymptotic behavior of the spectral density at 
low frequency is connected to the large time behavior of the covariance function. 
For the covariance function Ca,[3{T) which satisfles 

(3.24) C„,^(r)~L(r)|iT|i-"'3, ||t|1 ^ ex., 

where L{t) is a slowly varying function for large ||r||, i.e. L{ct) / L{t) — + 1 as ||t|| -^ 
oo for all positive constant c, Hardy-Littlewood-Karamata-Tauberian theorem (see 
e.g. [35j j36j) implies that the spectral density Sa^pioj) has the following asymptotic 
behavior 

So.A^) ^ c„,„^||c^|r'='-"L [^ , M ^ 0+, if aP < n, 

1 ^ [-^- 



where 



^n,a/3 



2"f^Tri 



(*) 



4. Gaussian Sheet with Generalized Gauchy Covariance 



We can also introduce another type of n-dimensional process X*a{t), called 
Gaussian sheet with generalized Cauchy covariance (GSGCC). 

Definition 4.1. The GSGCC is a centered Gaussian random field X* Jt) indexed 
by two multidimensional parameters a = (ai, . . . , a„) and (3 = {Pi, . . . , /9„), with 
ai € (0, 2], Pi > for alH = 1, . . . , n, and with covariance given by 

(4-1) C*,ir)^{x*^{t + r)X*At)) 

n n 

4=1 1=1 

where Ca^ i3i{Ti) is the covariance of one-dimensional GFGCC indexed by ai and 

P^■ 

Heuristically, one can regard X^ g{t) as the product of n independent one- 
dimensional GFGCC processes Xai^pi{ti). However, since the product of inde- 
pendent normal random variables is not a normal random variable, we shall not 
write X* Jt) — n"=i Xai,i3i{U)- Nevertheless, it is true that for any integer to > 1, 



■II 
E{[xtAi)Y)^X{E{[X^^^pMr). 



1=1 

The GSGCC ^^ ^(O i^ ^'^ anisotropic Gaussian random field and it does not 
satisfy the definition of local self-similarity given in Definition 12.21 However, we 
can show that it is lass and has tangent field according to Definition 12.31 Assume 
that a\ — . . . = a„i^ < a„i„+i < . ■ . < «„. Namely, for some to.q, G {1, . • . ,n}, 
ai = ^2 = ••■ = ctrjj.Q = mina, and for all i > ma -I- 1, a^ > mina. Here 
mina = min{ai, . . . , q;„}. 
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Proposition 4.2. GSGCC is lass with tangent field Ta^f3{t), which is a stationary 
centered Gaussian field with covariance 

rna 

1=1 
Proof. From definition, it is easy to see that 

n 

(4.2) C*p{r)^l[{l-(3,\nr[l + 0\r,n} as |r.| -- 
1=1 

Therefore, 

(4-3) C#^(r) = l-E/5^l^^l™" + E l^^r+E 

i^l \'i— 771^ + 1 i—1 

Consequently, 

' A,^X*^{t) A,,X*^{t)\ 



+ 



^ 1 2 min a 
Ti\ 



lim 



r.+ \ pniina/2 frmina/2 



: lim 






D 



From this proposition we see that in general, the tangent field Ta^p{u) of the 
GSGCC X* Jt) defined by Definition 12 . 31 are uncorrelated in some of the directions 
of M". It doesn't fully capture the locally self-similar property of X* Jt). There 
is another measure of local self-similarity, called local multiple self-similarity that 
better describe this kind of process. Recall that 19J a stochastic process X{t) is 
called multi-self-similar (mss) of index H E M" if and only if for any c G M" , 

X (Citi, . . . ,C„i„) =d Cj V . .C^"X{ti, . . . ,tn). 

It is obvious that a H-mss field is a [X)r=i -^d'SS field, but in general a self-similar 
field is not multi-self-similar. The notion of Iss can be generalized accordingly. 

Definition 4.3. A centered stationary Gaussian field is n-ple locally self-similar 
of order a/2 if its covariance function C{t) satisfies for ||t|| -^ 0+, 



C{r) = Y[{A-B,\nr[l + 0{\nf% 



i=l 



for some Ai, Bi, Si > 0, 1 < i < n. 

In view of (|4.2p . it is easy to see that X* Jt) is n-ple locally self-similar with 
Ai = 1, Bi = Pi and 5i = a^. We also note that the local multi-self-similarity is 
equivalent to 

(4.4) X*p(t + cne,) - X*p{t) =d c"' [xt^t + r.e,) - X#^(t)] as n ^ 0+, 

for any c £ M+. Here e^ is the unit vector in the i,; direction. This can also be 
rephrased as n-ple lass. 
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Proposition 4.4. XJ^ n{t) is n-ple lass. More precisely, for every i = 1, . . . ,n, 

( X*{t + eu,ei)-X*(t)] ^ 

(4.5) ^lim ^^ '^^ ^^^ = ^pB^^i,{u,\ 

where -Bc(/2(u) is the one dimensional fractional Brownian motion of index a/2. 



The proof is the same as Proposition 12.41 

One should observe that the n-ple lass here is different from the lass defined in 
Definition 12.31 We take the limit process in each of the ii-direction separately and 



the limit process can be regarded as the partial derivative process of X* Jt) in 
the ij-direction. A disadvantage of this definition is that the limit process of an 
n-dimensional field is a one-dimensional process. Therefore, we will define another 
limit process that will better reflect local multi-self-similarity, as were considered in 
[21] and [S]. Given a stochastic process X{t),t € M", we define its total increment 
by u e M" at the point i e M" as 

n„x(t)= Y. (-i)"~^^-^''^(i 

<5e{o,i}" \ 

When n ~ 1, this is the same as the increment, i.e. DuX{t) — A„X(i). When 
n = 2, we have 

□(ui,«2)^(il' ^2) = X{ti +Ui,t2+ U2) - X{ti + UiM) - X{ti,t2 + U2) + X{ti,t2), 

which is also known as the rectangular increment. For general n, given t,u di R", 
the set of points 




n 



U + J2S, 

I «=i ) (5e{o,i}" 

are the vertices of the hyperrectangle with t and t + u a.s one of the main diagonals. 
By giving the vertex t-\-u a weight -1-1, the other vertices of the hyperrectangle can 
be given a weight -1-1 and —1 alternatingly so that adjacent vertices has different 
weights. DuX{t) is then the weighted sum of the field X{t) at the vertices of the 
hyperrectangle. 

Recall that the fractional Brownian sheet -8^,2(^1* ^ ^"^ of index a/2 e (0, 1)" 
is a centered Gaussian process with covariance 



1 " 

Btf2it)B*/2i-'^)) = ^U^i^r- + \sr' -\u 



It is well known that B*,^{t) is a self-similar field of order X]"=i Q^j/2, and multi- 
self-similar of order a/2. However, unlike the Levy fractional Brownian field (|2.6p . 
the fractional Brownian sheet is not a process with stationary increment. Never- 
theless, it is a process with stationary total increment, i.e., 

{auB*^^{t), u e M"} =d {n„B# 2(5), u e R"} , yt,se R". 

Now returning to the local asymptotic multi-self-similarity of GSGCC, we can 
show the following: 
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Proposition 4.5. 




= d 



nv2A 



i=l 



B^/M- 



Here e.u = I]"=i ejWjgj. 



Proof. If we heuristically write X* Jt) as nr=i -^ai ,Pi (^* ) ' then heuristically the to- 
tal increment n„X*a(i) can be written as the product of increments YY^^i A^.Xq.^^^ (i^), 
and the result follows from Proposition l2.4l 
For a more rigorous proof, notice that 

<5e{0,l}" ^e{o,i}" \i^l / 

n 

1=1 

n 

The result follows. D 

From this proposition, we see that the total increment process can capture the 
locally multi-self-similar property of the GSGCC better. It also shows that locally, 
GSGCC behaves similarly as the fractional Brownian sheet. One would tend to use 
Proposition 14.21 and Proposition 12.61 to conclude that the Hausdorff dimension of 
the graph of GSGCC over a hyperrectangle is 

d„,a := n + 1 - - min{ai, . . . , a„}. 

However, Proposition 12.61 cannot be applied here since Proposition 14.21 does not 
imply that the fractal index of GSGCC is min a. For the fractional Brownian sheet 
B*,^, Kamont [23 showed that the Hausdorff dimension of its graph is bounded 
above by dn,a- In [1], Ayache used wavelet method to show that the Hausdorff 
dimension of the graph of the fractional Brownian sheet B^.^ is indeed equal to 
dn,a- This fact was proved again by Ayache and Xiao [6] as a special case of a 
more general result on Hausdorff dimension of fractional Brownian sheets from R" 
to M.'^, using a different method. Therefore, it is natural for us to conjecture that 
the Hausdorff dimension of the graph of GSGCC is also dn^a- In fact, (|4.3|) implies 
that 

4*,M = ([<;3(< + ^)-</3W]') -0(||r||-'-) as ||r|| ^0+. 

By a well-known theorem (see e.g. [1]), this implies that the Hausdorff dimension 
of the graph of GSGCC is bounded above by dn,a- On the other hand, it is easy 
to show that 
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Lemma 4.6. Let C = n"=i['^»i ^»] ^^ '^ hyperrectangle in . 
ci and C2 smc/i </iat 



There exist constants 



x*.pit)~xipis)] )<c2j2\u-s,r, 



(4.6) ci^|i,-s.r'< 

2 = 1 

/or all t, s £ C. 

Then by adapting the proof in 52 for the lower bound on the HausdorfF dimen- 
sion of general anisotropic Gaussian fields with stationary increments, it can be 
verified that the graph of GSGCC over a hyperrectangle has Hausdorff dimension 
equal to dn,a- 

The condition for LRD can be generalized to random sheet and it becomes 



(4.7) 
Since 



/ K.W 



d"T = 



Cf/3M=n^".A(r.), 

i=l 

and Proposition 12.81 for n = 1 gives 

\Ca,,f3iiTi)\ dn = OO <F=^ < aiPi < 1, 



/ 

K+ 



we have 



Proposition 4.7. GSGCC is LRD if and only if for some 1 < i < n, < aipi < 1; 

it is SRD if and only if ai(3i > 1 for all 1 < i < n. 

In order to determine the asymptotic behavior for the spectral density of the 
generalized Cauchy sheet, we observe that the spectral density of GSGCC, which 
we denote by S^j^ Juj) can be expressed as products of spectral densities of one- 
dimensional isotropic GFGCC considered in section [31 Namely, 

n 

(4.8) Slp{u)^\{So,a3M)- 

4=1 

Therefore, the high frequency and low frequency behaviors of the spectral density 
S* Jlo) can be obtained from Proposition 13.21 and Proposition 13.31 respectively. 

Proposition 4.8. The high frequency limit of the spectral density »5*g(w) is given 
by 



St0i^)-ll^^{a^^f3^ 



OO, 



where 






z/a, e(0,2); 



TLi{ai,(ii;uji 



ift-i 



2''-r(A) 



^-I'^.i 



if aj 
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Proposition 4.9. The low frequency limit of the spectral density S* Jto) is given 
by 

n 

1=1 
where 

A(a,,A;w,)= sm — -— wr*'^' \ if aiPi < 1; 

TT Z 

C^{ai,(3i;uji) =- \ In- — - - /J^ (^'(A) + 7) -7^ , «/ "ift = 1; 

TT I \uJi\ J 

ir(i-)r(A^i-) 

/:,(ai,/3i;wi) = — — , if ail3i>l. 

Here we have used the fact that T{2z) = 2'^''-'^Tr--^T{z)T {z + i) ([23], no.l of 
8.335), V{z)T{l~z) = 7r/sm(7rz) and ^(1/2) = -7-21n2 (i23j, no.2 of 8.366). 

5. Lamperti Transformation of GFGCC and GSGCC 

In his seminal paper |33j , Lamperti introduced a transformation which provides a 
one to one correspondence between a self-similar process and a stationary process. 
For a stationary process X{t), t G E, we let 

(5.1) Y{t)^t"X{lnt), 

for t G M+, H > 0, and Y{0) = be its 7J-Lamperti transfrom. Then Y{t) is an 
iJ-self-similar (H-ss) process. Conversely, if {Y{t),t > 0} is H-ss with zero mean, 
then the inverse Lamperti transformation of Y{t) defined by 

(5.2) X{t) = e-"'Y{e'), t e R, 

is a stationary process. 

There are two ways to define extensions of Lamperti transformation to M" linking 
stationary random field to a self-similar random field. Given a stationary random 
field X{t),t G M", the first way to define its Lamperti transformation is, given 
H €R+, defined by 

(5.3) Y{t)^\\t\\"X{lnti,...,lntn), 
for t e M" . It is easy to show that 

Proposition 5.1. Let Y{t),t g K^ be the H -Lamperti transform of a stationary 
field X{t),t e M" defined by dO]). Then Y(t) is a H-ss field. 

Proof. For any c G M+ , we have 

Y{ct)=d\\ct\\"X{ln[cti],...M[ctn]) 

=d c^||t||^X(lnti + Inc, . . . , lnt„ + Inc) 

=dc"\\t\\"X{lnti,...,lntn) 

=dc"Yit), 

i.e. Y{t) is H-ss. D 
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The inverse to this Lamperti transformation transforms a random field Y(t) to 

(5.4) X{t) ^ (e2*i + . . . + e2*")-^/2r(e*S . . . , e*") 

for t eW\ For a field X(t), i e R" to be stationary, it must satisfy X{t+u) ^d X{t) 
for all u G M". This is an 7i-parameter family of conditions. However, in proving 
that Y{t) is i?-ss, we only use one parameter family of these conditions, i.e., we 
only look at those u of the form ui = . . . = u„. Therefore we cannot expect that 
the inverse Lamperti transform ()5.4p of a H-ss field is a stationary field. As an 
example, consider the Levy fractional Brownian field with index H, Bnit), which 
is a H-ss field. Define its inverse 7J-Lamperti transform Z{t) by (|5.4p . We find 
that its covariance is 

(Z(t + T)Z{t)) =i (e2(*^+-i) + . . . + e2(*"+-")) '"^^ (e^'^ + . . . + e^*")"^/' x 

- ([e*^+^^ - e*^]2 + . . . + [e*"+-" - e*-]^]. 

It is easy to verify that for n > 2, this expression is not independent of t G M". 
Therefore, Z(t) is not a stationary process. 

To make full use of the n-parameter families of symmetry in a stationary process, 
there is another definition of Lamperti transformation introduced by [19| . Given a 
stationary process X{t),t S M", and a multi-index H S M" , the second iJ-Lamperti 
transform of X{t), denoted by Y(i), is defined as 

(5.5) Y{t) ^ if 1 . . . tf"X(lnti, . . . ,lnt„) 

for t e M" . It was shown in [19 that 

Proposition 5.2. Let Y(t),t g M" be the H-Lamperti transform of a stationary 
field X{t), t e M" defined by ^^. Then Y(i) is a H-mss field. 

The inverse of this second Lamperti transformation transform a field Y(t), t e M" 
to X{t),teW, where 

(5.6) X{t) = e-2:r=i *.ff.Y(e*S . . . , e*"). 

It has a nice property, namely, as was shown in [19j : 

Proposition 5.3. Let X(t),t G M" be the inverse second H-Lamperti transform 
of a H-mss field Y(i),t € M" defined by ()5.6|) . T/ien X(t) is a stationary field. 

As a side remark, it can be shown that (j5.5p is essentially the unique (up to 
some multiplicative constants) transformation that takes a stationary process to a 
H-mss. However, if for c e M" , /3 G K^, fc.p{t),t G R" is a function of the form 

n 
fc,p{t)=Y.c4, 
i=\ 

then if c[l], . . . , c\m\ S R" , f3i, . . . ,f}m S M"*" and 71, . . . , 7^ S R"*" are such that 
/?i7i -|- . . . -|- Pmlm = ii/^, it is easy to verify that the transform 



X{t) ^ I]/eH,ft(t)^*X(ln<i, . . . ,lni„) 
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takes a stationary process to a H-ss process, but its inverse in general does not 
take a H-ss process to a stationary process. The first Lamperti transformation 
(|5.3p corresponds to m = 1, c[l]i = . . . = c[l]„ — l^f3i — 2 and 71 = H/2. 

Returning to the GFGCC Xa,fj{t) and the GSGCC X*p{t), Proposition O 

shows that their first Lamperti transforms, Yaji{t) and y^ait), are H~ss fields 
with covariances 

(Yc,At)Yo.As)) ^\t\\"\\s\\"{Xo.A^nti, . . . ,lni„)X„,;3(lnsi, . . . ,lns„)J) 



and 



H\\„\\H 



\tru 



1+ ^(Int, -Ins,)' 



a/2' 



^i=l 



y*AmtAs)) =m"M"Xl [i + {^^u - in^o 



Q,l-ft 



respectively. On the other hand. Proposition 15. 21 shows that their second Lamperti 
transforms, Yaji{t) and Yj^(t), are H-vuss fields with covariances 

-/3 



Y„,^wY„,^(s))=n^f'n^?' 



4=1 i=l 



i/2' 



1+ ^(Ini, -Ins, 



and 



"^t.mtA')) -Wtf'^f' [1 + ^^^t, - lnsO"T 



respectively. 

Next we consider whether the LRD property of GFGCC and GSGCC is preserved 
under the Lamperti transformations. First we note that the LRD condition (|2.9p 
can be extended to non-stationary field in the following way. 



Definition 5.4. Let 

(5.7) Rx{t,t + T) 



C{t,t + T) 



y/Cit + T,t + T)C{t,t) 



be the correlation function of a non-stationary Gaussian field X{t). Then the 
condition of LRD for the non-stationary Gaussian field X{t) is given by 



(5.8) 



/ \Rx{t,t + T)\d"T ^00. 



Proposition 5.5. 

A. The H-ss Gaussian field Yajj{t) is LRD for all a G (0, 2] and /3 > 0. 

B. The H-ss Gaussian field Y^g{t) is LRD for all a,f3 E M" satisfying Oi S (0,2] 
and Pi > 0, 1 < i < n. 

C. The H-mss Gaussian field Ya.pit) is LRD for all a G (0, 2] and /3 > 0. 

D. The H-mss Gaussian field Y*^(i) is LRD for all a, (3 E K" satisfying ai G (0, 2] 
and Pi > 0, 1 < i < n. 
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Proof. It is easy to verify that both the processes Ya^f^it) and Ya,p{t) have the same 
correlation function 



(5.9) 



RY^At + r,t) = RY^Jt + T,t) 



1+ E^ 



u 



a/2' 



T e 



^+' 



which is independent of H. In order to show that these processes are LRD, we have 
by condition (15. 8p . 

-/3 



|i?y^^(i + r,i)|d"r 



1+ E >" 



-f 



a/2 



Wu 



Wu 



l+^(ln[l + z.,]r 



a/2 



-0 



n \ "/2" 


-p 


n 


>:-? 




He- 


i=i / 




.1=1 



d"t;, 



where we have used the substitutions Ui — Ti/ti and Vi — ln(l + Ui). Clearly, when 
||w|| -^ oo, the integrand also approaches oo. Therefore, the last integral diverges for 
all a,/3, which is the condition for Ya^pit) and Ya.p{t) to be LRD. The statements 
for Y^Jt) and Y* At) are proved analogously. D 

Note that the LRD property of GFGCC and GSGCC is preserved under the 
Lamperti transformations, but the SRD property is not preserved. Here we have 
an example that the application of Lamperti transformation to a LRD stationary 
process (in this case the GFGCC and GSGCC) gives a (multi)-self-siniilar process 
with LRD. Examples of Lamperti transformation encountered so far relate either 
two short memory processes (for example, Ornstein-Uhlenbeck process and Brow- 
nian motion), or between a SRD process and a LRD process (in the case of fBm 
and its inverse Lamperti transformed process). For examples, one can show that 
the inverse Lamperti transformation of fractional Brownian sheet with LRD prop- 
erty gives rise to a stationary random sheet with short range dependence, and the 
stationary field associated with the Levy fractional Brownian field is a stationary 
field with SRD. 

Recall that Levy fractional Brownian field is essentially the only self-similar 
Gaussian field with stationary increments |47j . Similarly, one can show that frac- 
tional Brownian sheet is essentially the only multi-self-similar Gaussian random 
field that has stationary total increments. Hence Va,/3(i), ^a /3(0' ^a,p{t) and 
Y* Jt)^ the self-similar and multi-self-similar fields associated with GFGCC and 
GSGCC, do not have stationary increments or total increments. There exists a 
weaker stationary property known as asymptotically locally stationarity, which re- 
quires the field to be stationary in the limit ||r|| — > 0"*". 

Definition 5.6. A centered Gaussian random field X{t) is said to have asymptot- 
ically locally stationary increment if and only if as ||t|| — > 0+, the variance of its 
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increment (J^{t) — ( [ArX{t)] \ is independent of t. More precisely, 

(5.10) a^{T)=f{T)+g{t,T), 

where /(r) is independent oft, and for any fixed i, g(t, r) = o{f{T)) as functions of 
T. Similarly, a centered Gaussian random field X{t) is said to have asymptotically 
locally stationary total increment if and only if as ||t|| -^ 0+, the variance of its 

total increment (T^(t) — / [n^X(t)] \ is independent of t. 



For the random field yQ./3(t), we have 

2\ "/2' 



= \\t + T 



]2H 



\tf"-2\\t + T\\"m" 



1+ E 1- 






Using 



\\t + r\\ = Vih + ny + ... + (t„ + T„)2 =||t|| (^1 + ^ ^g^/^""^ + o(lklP)) 



1/2 



ill 1 + 



and 



h?] 


.2\"/2" 


-/:! 



l + ^^ + 0(||r|p) 



a/2' 



-ft 



2 \ "/2' 

+ 0(||rf) 



-/3 






\i=l 



i/2 



+ 0(||r|r+^) 



q/2 



O IIt 



I min{2a,Q:+l} 



we find that as ||t|| -^ 0^, 



(5.11) ( [ArY^At)?) =\\t\r ( 1 + 2H^^^) + \\t\r 



-2||t||2» [l + H 






iitip 

II^P 

min{2,2a,a+l}^ 

2\ "/2 



1-ME 



^i=l 



a/2N 



^4=1 



(^ n|^||min{2,2a,Q+l} 
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It is obvious that the leading term is not independent of t for any H and a, unless 
when n ~ I and a = 2H. Therefore for n > 2, Yajj{t) does not have asymptotically 
locally stationary increments. Similarly, for the fields Y^jjit), Y^bW ^^^ ^a Bi^)^ 
one can show similarly that as |1t|1 -^ 0+, 

2\ "/^ 



(5.12) {[ArY^Atf) =2/3tr^ . . .e- E 



Q n|^||min{2,2a,a+l} 



(5.13) 



^rY*^it)] )=2iJtf^E'3^ F +o(||^|jmi„{2,2«.,«. + l}^^ 



(5.14) /[A,Y#^(t)] U2i?^^..e"E^^ F +0(||T|r-^2-2"— +i>). 

Therefore for n > 2, none of the fields Ya.fj{t), Y*Jt) and Y* Jt) have asymp- 
totically locally stationary increments. When n = 1, all the fields Ya^/3{t), Yc,/3(<), 
y^o(t) and Y* p{t) are actually the same, and they have asymptotically locally sta- 
tionary increment if and only if a = 2_ff , in which case the variance of the increment 
cr|(T) behaves like 

aUr) ~ 2/3|rr \r\ ^ 0. 

Next wc consider the total increments. Since increment is the same as total incre- 
ment when n — I, we only need to consider n >2. Using similar computations as 
given above, one can verify that for n > 2, the fields Ya^i3{t), Ya.jsit) and Y^Jt) do 

not have asymptotically locally stationary increments, but Y* Jt) have if a = 2H. 
We show the computation of the latter case here. By definition, 

2\ 



^rYl,{t) 



E E (-1)^^^^'+^^^- 

<5e{o,i}" r,e{o,i}" 

( </3 ( * + E ^^^A ^t(3 ( ^ + E ^^^^' 



i=i (5,,J7i)e{04}2 

ai\ -Pi ~\ 



=n«f--2(«. +'-.)"•«.''' 1 






=n 2M"' 



o 



In 

niin{2,2ai,Q:i + l} 



+ {U+Ti) 



2 Hi 



-OM|t||^-i°'+'^ 
where 6 = min{2 — ai,ai, l}"^i- If a = 2_ff, then as ||tJ| -^ 0+, 

(^[n,Y*^{t)]''j ^2-f[mnn- 

Therefore, we have verified the following: 



4 = 1 
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Proposition 5.7. The total increments of Y* Jt) are asymptotically locally sta- 
tionary if a = 2H . 

Next we consider the tangent fields (Definition 12. 3p of the Lamperti transforms 
of GFGCC Ya.f}{t) and Yc,./3(i). We have 

Proposition 5.8. 

A. The field Ya.p{t) is lass of order a/2, with tangent field at t ^ M" being 



y2:0iitrs„/,(^,...,^) 



B. The field Ya^f3{t) is lass of order a/2, with tangent field at t ^ M" being 



.i=l 



Hi 



-Da/2 I ^, • • • , ^ I • 



Proof. Using the formula 

{^e^Xa..fi{t)l^evYa,.f}{t)) 

we obtain immediately from (|5.1ip that as e — > 0, 

'Ae„ra,/3(i) Ae«r„,^(i)' 




E 



ti 



q/2' 



Notice that this is up to the factor 2/3||t|p^, the covariance of the scaled Levy 
Brownian field 

The statement for Yo.,ii{t) is proved similarly. 



u e 



D 



In fact, one can also deduce from (|5.1ip and (|5.12p that the local fractal index of 
the fields Ya^i3{t) and YQ,^^(i) are both equal to a/2, and the Hausdorff dimension 
of their graphs over a hyperrectangle is n + 1 — a/2. 

For the tangent fields of the Lamperti transforms of GSGCC Y*Jt) and Y* „{t), 
we recall that given a S (0,2]", mina — minjai}"^]^, and we assume WLOG that 



mm a — ai 



< a. 



rria + l 



< 



< a„ for some 1 < rrin < n. Then as in 



Proposition l5.81 we can use (|5.13p and (|5.14p to show that 

Proposition 5.9. 

A. The field Y*p{t) is lass of order mina/2, with tangent field at t £ R" being 



\t\\"T, 



a, 13 



Wl 






30 



S.C. LIM^ AND L.P. TEO^ 



where the field Taj3{u)^ u £ R" is defined in Proposition \4.S\ 

B. The field Y*g(i) is lass of order mma/2, with tangent field at t £ R" being 



.1=1 



Hi 



T. 



a, 13 









From Propositions 12.41 and 15.81 and Propositions 14.21 and 15.91 we find that the 
tangent fields of GFGCC and GSGCC are related to the tangent fields of their 
Lamperti transforms by some change of variable formulas. Namely, if X{t) has 
tangent field T{u) at i G R", then its first Lamperti transform Y[s) has tangent 
field 

ll-sirrf^,...,^ 

\Si Sn 

at s e R" , and its second Lamperti transform Y(s) has tangent field 



n 



.Hi 



T 



Ui 


u 


Si 


• 1 

s, 



at s E R" . We also notice that the order of self-similarity of the fields ^a./slO ^-nd 
y^ait), H, is in general different from their order of local asymptotic self-similarity, 
being a/2 and mina/2 respectively. 

Since the stochastic process Y^^(i) resembles the GSGCC Xfjt) in the sense 
that both their covariances can be written as products of covariance of one-dimensional 
processes, we can consider the limit of the total increments oiY* p{t). It is easy to 
obtain as in Proposition 14. 5i using the n = 1 case of Proposition 15. i 



that 



Proposition 5.10. 



Here e.u 




— d 



,H,-ia,/2) 



.4=1 



sf/^H- 



6. Concluding Remarks 

We have studied some of the basic properties of GFGCC and GSGCC, and their 
associated Lamperti transforms. The asymptotic properties of the spectral densities 
of GFGCC and GSGCC are considered. One expects the separate characterization 
of fractal dimension and long range dependence for GFGCC and GSGCC will pro- 
vide more flexibility in their applications to modeling various surfaces and images. 
In the one-dimensional case, GFGCC has been applied to model the Havriliak- 
Negami relaxation law [38] . The estimations of parameters for stationary Gaussian 
processes have been widely studied. Some of these estimations can be adapted for 
GFGCC and GSGCC, and they are crucial to the applications. Further applications 
in spatial-temporal processes are possible if GFGCC is extended and modified to a 
space-time field to include non-stationarity and anisotropy [4T1 [32l [9l [20l l48l [211 [T3] . 
We hope to apply results obtained in this paper to model various physical systems 
such as thin film surfaces in semiconductors [43l |46] , surface ocean waves [^ , geo- 
logical morphology ^3\, etc, in a future work. 
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